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Evolution of Supernova Remnants Expanding out of the Dense 
Circumstellar Matter into the Rarefied Interstellar Medium 
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Abstract 

We carry out 3D-hydrodynamical calculations for the interaction of expanding supernova ejecta with 
the dense circumstellar matter (CSM) and the rarefied interstellar medium (ISM) outside. The CSM 
is composed of the stellar wind matter from the progenitor in its pre-supernova phase, and assumed to 
be axially symmetric: more matter around the equator than in the polar direction driven by rotation 
of the progenitor. Because of high density of the CSM, the ionization state of the shock-heated ejecta 
quickly becomes equilibrium with the electron temperature. When the blast wave breaks out of the CSM 
into the rarefied ISM, the shocked ejecta cools rapidly due to adiabatic expansion, and hence an over- 
ionized/recombining plasma would be left. The ejecta is reheated by the second reverse shock due to the 
interaction with the ISM. We calculate the emission measure of the supernova remnant (SNR) along the line 
of sight, and find that the over-ionized plasma appears to be bar-like with wings in the edge-on (equatorial 
view), while shell- like in the face-on (polar view) geometry with respect to the rotation axis. The hot 
gas heated by the blast wave exists in the outermost region of the SNR with a nearly complete shell, but 
the X-rays therefrom are too faint to be observable. Thus, depending on the viewing angle, the SNR of 
the over-ionized plasma would exhibit center-filled morphology in X-rays, like W49B, a mixed-morphology 
SNR. The bar-like structure is swept out by the second reverse shock and disappears eventually, and then 
the SNR becomes shell-like in both the equatorial and polar views in the later phase of the evolution. 

Key words: hydrodynamics - stars: circumstellar matter - ISM: supernova remnants - methods: nu- 
merical 



1. Introduction 

Evolution of a supernova remnant (SNR) with no neu- 
tron star in its center is described essentially by the shock 
waves propagating outward (blast wave) into the ambi- 
ent matter and inward (reverse-shock wave) into the su- 
pernova ejecta. In young SNRs, the shocked matter is 
heated up to temperatures higher than 10 7 K and forms 
an optically-thin hot plasma. This plasma emits X-rays 
with many emission lines of highly- ionized heavy elements. 
Since the time-scale of ionization by electron impact is 
longer than that of shock-heating of electrons, X-ray spec- 
tra of young SNRs exhibit ionization states lower than 
those expected from the electron temperatures, which is 
called under-ionized or ionizing. The ionization state is 
represented by ionization temperature T z , the plasma tem- 
perature in collisional ionization equilibrium having the 
relevant ionization state. Thus, many of young SNRs show 
T z < T c (under-ionized/ionizing) for the electron temper- 
ature T c . 

On the contrary to the above standard scenario, 
Kawasaki et al. (2002) found an over-ionized plasma (T z > 
T c ) in the ASCA observation from a mixed-morphology 
SNR, IC443. They proposed that thermal conduction 
from the hot interior of the remnant to the cold exte- 
rior can explain the over-ionized plasma. Kawasaki et al. 



(2005) further investigated other five mixed-morphology 
SNRs, W49B, W28, W44, 3C 391 and Kes 27, and found 
that W49B shows T z > T c as well. 

Suzaku detected clear spectral structures of radiative- 
recombination continua, direct evidence for recombining 
plasmas, in the X-ray spectra of IC443 (Yamaguchi et al. 
2009) and W49B (Ozawa et al. 2009). Radiative recom- 
bination continua, and hence an over-ionized plasma is 
also found in another mixed-morphology SNR, G359.1-0.5 
(Ohnishi et al. 2011). Although the samples are still lim- 
ited, the over-ionized plasmas are all found in the mixed- 
morphology SNRs associated in or near the star-forming 
complex with HII regions and molecular clouds. Hence 
these SNRs are likely due to core-collapse supernovae; the 
over-ionized plasma is possibly related to the massive pro- 
genitor and/or its environments. 

One possibility is that the ionization by photons of a few 
tens keV at the initial phase of explosion of massive stars, 
e.g., possible X-ray flash or afterglow of 7-ray bursts, is 
responsible for the over-ionized state. Another possibility 
is rarefaction in the adiabatic evolution phase of the SNR; 
when the shock wave breaks out of the dense circumstel- 
lar matter (CSM), like the progenitor's stellar wind, into 
the rarefied interstellar medium (ISM), the shock- heated 
electrons would rapidly cool due to adiabatic expansion, 
leaving highly ionized states (Itoh & Masai 1989). 
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In the present paper, we investigate the latter possibil- 
ity. Itoh & Masai (1989) calculated SNR evolution and 
radiation for a spherically symmetric structure. We ex- 
tend this work to non-symmetric structures as are found 
in mixed-morphology SNRs. Unlike Itoh & Masai (1989), 
we do not take ionization or radiation processes into ac- 
count; we assume the relevant SNR evolution is fully adi- 
abatic. We describe the model in the following section, 
calculation and results in Section 3, and discuss the re- 
sults in Section 4. We give a summary in section 5 and 
notes about electron temperature in Appendix. 

2. Model 

2.1. Circumstellar matter 

In a core-collapse supernova, the massive progenitor 
(red or blue supergiant or Wolf-Rayet star) is expected 
to blow a strong stellar wind in its pre-supernova phase 
and forms dense CSM around the progenitor. If the stel- 
lar wind is spherically symmetric with a mass-loss rate M 
and wind velocity v w , the CSM density p at a distance r 
from the progenitor is given by 

p _ M 
4irv w r 2 

For a while before supernova explosion, the stellar wind 
ceases. Then the stellar wind forms a thick shell. The 
inner and outer radii of the shell are given by i?; n = v w t e 
and -R out = v w (t w +t c ), respectively (see Figure 1), where 
time t w and t e are the duration of the wind activity and 
the elapsed time after the cease of the wind, respectively. 

In realistic stellar winds, more matter around the equa- 
torial plane than the polar direction may be accumulated, 
because of the rotation of the progenitor; disk-like winds 
of Be stars for an extreme example. We thus modify Eq. 
(1) as, 

M a 



4nv w x 2 + y 2 + a 2 z 2 ' 
Rin<{x 2 + y 2 + a 2 z 2 )^ 2 <R 



(2) 



with an anisotropy parameter a, which is the ratio of the 
major to minor axis of an equi-density surface of the wind 
matter. Hereafter, R ln and i? ut are referred to the radii 
on the equatorial plane. This distribution gives a 2 times 
higher density on the equatorial plane than in the po- 
lar direction at the same distance. For the progenitor of 
SN 1987A, Blondin & Lundqvist (1993) suggested that 
the ratio of equatorial to polar mass-loss rate was at least 
20 during its red supergiant stage. This value corresponds 
to a ~ 4.5. In our calculations we adopt the values from 
1 to 3. 

We calculate the evolution for the following five models: 

Al Isotropic CSM same as Itoh & Masai (1989) for a 
reference. 

A2 Anisotropic CSM concentrated around the equato- 
rial plane. 

A3 Same as A2, but more mass around the equatorial 
plane. 



Bl Same as A2, but with a lower density ISM. 
B2 Same as Bl, but with a higher wind velocity. 

The numerical values for each model are shown in Tabic 
1, where values of i?j n , i? ut, Pin and p out are those on 
the equatorial plane. The mass-loss rate is assumed to be 
5 x Kr 5 Af yr -1 , the same as Itoh & Masai (1989), for 
all the models. 
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Fig. 1. Schematic picture of initial configurations for the in- 
teraction of the supernova ejecta with CSM. The broken and 
solid lines of the ejecta represent the radii of the core and 
envelope, respectively, and the outer two solid lines represent 
the inner and outer boundaries of the CSM to the ISM. 



2.2. Supernova ejecta 

We focus on the effect of the CSM, particularly for 
anisotropic CSM, on the early phase evolution of SNRs. 
For the supernova ejecta, we simply assume spherical dis- 
tribution with a core of radius R c and extended envelope 
of radius R e 



3M 



j, as (Truclove & McKee 1999) 
l-n/3 



iirv^t 3 1 - (n/3) 



x 



(r/(v ei t)Y 



U'r 



for r < R c 
for R c <r< Re 



(3) 



where 



(2E e 



1/2 



1/2 



,-(3-n) 



■n/3 



1/2 



w, 



-(5-n) 



■n/5. 



Wr 



is the expansion velocity at i? e j- Afej and E e j are the 
ejecta mass and the explosion energy (kinetic energy), re- 
spectively. These are assumed to be E C j = 2x 10 51 erg and 
M e j = 10 Mq. The other parameters n and w c = R c / R e] are 
taken to be 6 and 0.49, respectively. All these values are 



the same as Itoh & Masai (1989), and hence v e 
cm s _1 . 



:8.5 x10 s 



3. Calculation and Results 

We numerically solve the Euler equations utilizing the 
athena3d code (Stone et al. 2008): 
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Table 1. Parameters of circumstellar matter and interstellar medium 











CSM 






ISM 


Model 


a 


i? in (10 ib cm) 


R ont (10 iO cm) 


Pin (amu cm -3 ) 


Pout (amu cm - - 5 ) 


v w (cm s x ) 


piSM (amu cm" J ) 


Al 


1 


2.0 


92.6 


3.4xlO b 


1.6x10* 


10 b 


1.6X10" 1 


A2 


2 


2.0 


92.6 


6.9xl0 5 


3.2xl0 2 


10 6 


1.6X10" 1 


A3 


3 


2.0 


92.6 


l.OxlO 6 


4.8xl0 2 


10 6 


1.6X10" 1 


Bl 


2 


2.0 


92.6 


6.9xl0 5 


3.2xl0 2 


10 6 


1.6xl0~ 2 


B2 


2 


20 


926 


6.9xl0 2 


3.2X10- 1 


10 7 


1.6xl0~ 2 



Mass loss rate M is 5 x 10 5 M Q yr 1 in the all models, and Pi„ = p(r — i?i n ), p out = p(r — R ou t). 



dv , 1 ._. 

— + (v-V)v = --Vp, (5) 

m- 2 +-)}+v.{p.(^ 2 +u +Q}=o, 

where p, p, v and u are the mass density, the pressure, 
the fluid velocity and the internal energy per unit mass, 
respectively. The equation of state is written as p = (7 — 
l)pu with 7 = 5/3, and the mean molecular weight is taken 
to be 0.5. 

In the hydrodynamical evolution, we obtained mean 
temperature T of ions (protons) and electrons weighted by 
their number density. If no plasma mode works directly on 
electrons at the shock front, the ion temperature 7] rises 
faster and then the energy of ions is transferred to elec- 
trons. We calculate 7] and T e from T, assuming the trans- 
port by Coulomb collisions (Masai 1994; see Appendix). 

3.1. Model Al 

Although we include neither ionization nor radiation 
and simply assume the adiabatic evolution, the evolution 
of model Al is in fairly good agreement with Itoh & Masai 
(1989). The blast shock and reverse shock heat the CSM 
to T > 10 8 K, and the ejecta to T > 10 7 K, respectively. 
The blast wave expands as oc i - 88 , while r D oc t - 87 in 
Itoh & Masai (1989), where r D is the blast-wave radius. 

The blast wave breaks out at 41 yr, 2 years earlier 
than in Itoh & Masai (1989) who take radiation loss into 
calculation. After the break-out, the blast wave rapidly 
expands adiabatically, and hence the temperature and 
density of the shocked matter decrease approximately as 
T oc t~ 2 and p oc t~ 3 for 7 = 5/3. 

When the pressure of the shocked CSM becomes below 
that of the shocked ISM, the second reverse shock occurs 
to propagate inward and reheats the CSM and ejecta. The 
second reverse shock reaches the ejecta at 830 yr, 120 years 
earlier than in Itoh & Masai (1989). 

3.2. Models A2 and A3 

The evolution is basically the same as model Al, but 
is dependent on direction: more matter on the equato- 
rial plane than in the polar direction. With increasing 
the value of a, the break-out occurs earlier in the po- 
lar and later in the equatorial direction than in model 
Al. The former is due to a shorter distance to the CSM, 



and the latter is higher density of the CSM than those in 
model Al. In addition, the second reverse shock reaches 
the ejecta earlier in the polar direction and later in the 
equatorial direction than in model Al. In model A2, the 
break-out occurs at 17 yr in the polar direction, and at 
46 yr in the equatorial direction as shown in Figure 2. The 
mean temperature T of the shocked ejecta turns to rise at 
340 yr and 1100 yr, in the polar and equatorial directions, 
respectively. 
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Fig. 2. Shocks in the equatorial direction of model A2 before 
and after the blast- wave break-out, as functions of the elapsed 
time after explosion. The upper, middle and lower solid lines 
represent the radii of the blast wave, second reverse shock and 
reverse shock, respectively. The upper and lower broken lines 
represent the outer radii of the CSM and ejecta, respectively. 

3.3. Model Bl 

The break-out occurs at the same time as model A2 
since the CSM distribution is the same as model A2. The 
difference of model Bl from model A2 is the lower density 
of the ISM by an order of magnitude. The velocity of the 
blast wave after the break-out is slightly higher than that 
in model A2 by a factor of ~ 1.4 in the equatorial direction. 
As a result, the pressure of the shocked ISM is about one- 
fifth of that in model A2, and the second reverse shock 
occurs later. The temperature T of the shocked ejecta 
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drops at 23 yr and 68 yr by the break-out, but turns to 
rise at 810 yr and 1800 yr in the polar and equatorial 
directions, respectively. 

34. Model B2 

The difference from model Bl is the lower density of 
the CSM by three orders of magnitude due to higher wind 
velocity by an order of magnitude. Even at such low den- 
sities, the break-out occurs at 170 yr and 450 yr in the 
polar and equatorial directions, respectively. Similarly to 
the other models, the temperature T of the shocked ejecta 
drops at 280 yr and 770 yr, but turns to rise by the sec- 
ond reverse shock at 920 yr and 2400 yr in the polar and 
equatorial directions, respectively 

4. Discussion 
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In Table 2 we summarize the temperatures and blast- 
wave radii at characteristic epochs, t^ and t^ (see below), 
for each model. 

Dynamical evolution 

Before the break-out, the evolution of SNR can be de- 
scribed by a self-similar solution, and the blast wave ex- 
pands as rb oc t a . Almost independently of anisotropy of 
the CSM (models Al-3), the values of a are close to each 
other both in the polar and equatorial directions, as seen 
in the upper panel of Figure 3. After the break-out, a 
values are slightly larger in larger a models, but approach 
the value of model Al (a = 1) with time. This implies 
that the blast wave approaches spherical symmetry with 
time, and the ratio of the blast- wave radius of the polar to 
equatorial direction becomes less than ~ 1.2 in 1000 years 
in any model of ours here. 

If the density difference across the interface between 
the CSM and ISM is large enough (models Al-3 and Bl), 
the blast wave gets faster just after the break-out by a 
factor of ~ 2, almost independent of the density difference. 
On the other hand, if the density difference is as small 
as model B2, the velocity increases by a factor of ~ 1.4. 
The expansion velocity of a spherical fluid initially at rest 
asymptotically approaches its maximum value of 2c s /(7 — 
1), where c s is a sound speed of the fluid (Zel'dovich & 
Raizcr 1966). Similarly, denoting the shock velocity by 
Vs, we may have the maximum expansion velocity after 
the break-out, as 

^k 2 + F s 2 (t b )] 1/2 ^2Ws(tb), (6) 
7-1 

where Vs(£t>) is the shock velocity of the blast wave im- 
mediately before the break-out at t = t^. Here we take 
the Rankine-Hugoniot relation of strong shocks into ac- 
count. This estimate is consistent with the result of our 
hydrodynamical calculation. The high shock velocity and 
rarefaction caused by the break-out may be in favor of 
particle acceleration, because the former boosts the maxi- 
mum energy of particles, and the latter works for particles 
to become non-thermal. 



Fig. 3. Radius and velocity of the blast wave in the polar 
direction for models Al (solid lines), A2 (dashed lines) and 
A3 (dotted lines), as functions of elapsed time after explosion. 

4.2. Temperature and Ionization state 

As demonstrated by Itoh & Masai (1989), rarefac- 
tion by the blast-wave break-out rapidly cools the shock- 
heated matter to be an over-ionized/recombining plasma 
of T c <T Z . This recombining plasma state lasts until the 
second reverse shock reheats the matter to > T z . Although 
ionization is not calculated simultaneously with hydrody- 
namics, we can analyze the ionization state, i.e., ioniz- 
ing, equilibrium or recombining, of the shocked matter 
from the electron temperature and density. The ioniza- 
tion state thus evaluated is consistent with Itoh & Masai 
(1989), considering the variation of epochs and duration 
that depend on the CSM/ISM models, as discussed in the 
previous sections. 

The upper panels of Figures 4a (model Bl) and 4b 
(model B2) show the evolution of the ion and electron tem- 
peratures averaged over the shocked matter of the ejecta. 
Before the break-out, T c becomes nearly equal to 7] due 
to high densities of the CSM. Also after the break-out and 
until arrival of the second reverse shock, T e is nearly equal 
to Ti, because cooling is due to adiabatic expansion. By 
the second reverse shock, Ti rises faster than T c , and then 
T e rises through the energy transport from ions, as shown 
in the upper panels of Figure 4. 

In the lower panels, we show the average density of the 
shocked ejecta by the solid line. Ionization by electron- 
impact becomes equilibrium with time-scale r given by 
pr ~ 10 12 amu cm -3 s, almost independent of the tem- 
perature or ion species (Masai 1994). We also plot this 
relation as p C rit = 10 12 t -1 amu cm" 3 with the dashed line 
in the lower panels. Well before the break-out, since the 
density is high enough (p > p cr it) due to the presence of 
CSM, ionization quickly reaches its equilibrium at T e , i.e. 

As the SNR expands, the average density of the shocked 
matter decreases approximately as poet" 2 . When rarefac- 
tion occurs by the break-out, p decreases faster as oc t -3 
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Table 2. Characteristic Epochs 







Break-out 




T 

± z 


— T e Decoupling 


Model 


*b (yr) 


T (10 Y K) 


i? b (pc) 


*d (yr) 


T (10'K) 


R h (pc) 


Al 


41 


1.3 


0.30 


180 


0.20 


1.7 


A2-e 


46 


1.2 


0.30 


210 


0.18 


1.8 


A2-p 


17 


1.6 


0.15 


80 


0.16 


0.99 


A3-e 


47 


1.4 


0.30 


220 


0.21 


1.8 


A3-p 


11 


0.8 


0.10 


60 


0.11 


0.81 


Bl-e 


46 


1.2 


0.30 


190 


0.23 


2.0 


Bl-p 


17 


1.7 


0.15 


80 


0.10 


1.1 


B2-c 


450 


1.2 


3.0 


70 


1.7 


0.64 


B2-p 


170 


1.1 


1.5 


36 


1.2 


0.38 



The characters "e" and "p" attached to Model A2-B2 mean the equa- 
torial and polar directions, respectively. 




Fig. 4. Model (a) Bl and (b) B2: Averaged temperature and density of the shocked ejecta in the equator (left) and polar (right) 
directions, as functions of the elapsed time after explosion. Upper: The solid and broken lines represent the electron and ion 
temperatures, respectively. Lower. The solid and dashed lines represent the density and p C rit = 10 12 t~ 1 amu cm -3 , respectively, 
and the thin-dotted vertical lines represent the decoupling epoch t d (see text). 



(see Section 3.1). Therefore, p becomes below p cr it oci" 1 , 
as seen in the lower panels of Figure 4, and T z is decoupled 
from T c . In other words, recombination no longer follows 
the rapid decrease of the electron temperature, and the 
ionization state freezes roughly at T z ~ T c (td), where id is 
the epoch at which p = p CT it- 

The decoupling epoch t<j is also shown in the lower pan- 
els of Figure 4 with the thin-dotted vertical line. In model 
Bl, ta ~ 190 yr in the equatorial direction. Note that 
id > tb and T c (id) < T e (tt>) in model Bl, where T e (ib) is 
the electron temperature immediately before the break- 
out. Models Al-3 show the similar behavior to model Bl 
described here. 

In model B2, unlike other models or Itoh & Masai 
(1989), the CSM is located away from the progenitor, and 
hence the density is low. Therefore, the break-out occurs 
later but p becomes below p cr ; t earlier than in model Bl. 
As seen in Figure 4b, id < % in model B2, and id ~ 70 yr 
in the equatorial direction. In fact, as is seen in the up- 
per panel, T e before the break-out is nearly constant at 
~ ?o(fb) in model B2. As a result, the ionization tem- 



perature of the recombining plasma is higher in model B2 
than that in model Bl. 

When the second reverse shock reaches the ejecta, the 
temperature turns to rise, as seen in the upper panels of 
Figure 4. If the electron temperature exceeds or becomes 
comparable to T z ~ T (id), the ionization state turns to be 
ionizing or nearly equilibrium. This is the case for models 
Al-3 with a higher ISM density. In model Bl, though the 
ISM density is low, T c exceeds T z in 3000 years owing to 
the dense CSM and lower T z . On the other hand, in model 
B2, the density is too low for the second reverse shock to 
raise T c > T z in such short time, and the over-ionization 
state remains much longer. If 7] becomes higher than 2T Z 
at the shock front, T c ^ T z may be attained by the energy 
transfer from ions. Since T; oc V s 2 (reverse shock) decreases 
with the age, however, the recombining state would last 
yet for thousands years. 

After decoupled from the electron temperature, the 
ionization temperature decreases by recombination with 
time-scale r. We note here that the recombination time- 
scale ~ 10 4 yr (Masai 1994) is longer than the age of over- 
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ionized SNRs: - 4000 yr for IC443 (Troja et al. 2008) and 
1000 - 4000 yr for W49B (Pye et al. 1984, Smith et al. 
1985, Hwang et al. 2000). 

4- 3. X-ray emission measure distribution 

All the SNRs with over-ionized/recombining plasma so 
far observed are mixed-morphology SNRs, which exhibit 
center-filled X-ray emission and shell-like radio emission. 
Hence, we investigate distribution of the X-ray emission 
measure for SNRs with anisotropic stellar wind. We inte- 
grate the density square along the line of sight, as / p 2 dl, 
and show the map (black lines) for model B2 in Figure 5. 
Here the shocked matter of T > 10 6 K, which is responsible 
for X-ray emission, is taken into calculations. 

At ages 980 yr and 1800 yr, we can see a bar-like struc- 
ture (horizontal) with diffuse wings in the equatorial view, 
while a thin shell in the polar view. The emission measure 
is dominated by the shock-heated ejecta, which is a recom- 
bining plasma as seen from Figure 4b. The temperature 
dependence of the line emissivity is weaker in recombin- 
ing than in equilibrium, and the emission-measure map 
reflects roughly the X-ray surface brightness. 

The bar structure is clearer at younger ages after the 
break-out, and diffuses out gradually with time. At 
1800 yr, the second reverse-shock front appears to sur- 
round the bar/wing structure, and thereby the bar-end 
structure is being distorted. The second reverse-shock 
propagates inward and eventually sweeps the bar/wing 
structure out. Hence, the late-phase SNR approaches 
shell-like both in the equatorial and polar views, as seen 
at 10000 yr in Figure 5, where a little elongated shape 
with a narrow middle part is still seen in the equatorial 
view. 

W49B shows a bar-like structure as well as a recom- 
bining plasma state (Ozawa et al. 2009). The bar and 
its east-end observed with Chandra (Keohane et al. 2007) 
look similar to our calculation shown in Figure 5. Thus we 
suggest that W49B is the case of nearly equatorial view 
for the supernova remnant of a massive progenitor which 
exploded in its past stellar wind. It should be noted that, 
in our calculation, the structure of the bar-end is formed 
by the reverse shock, not by the collision with a molecular 
cloud near the bar-end. 

The grey lines in Figure 5 represent the low-level emis- 
sion measure, much smaller than that represented by the 
black lines. One can see the blast-wave front in the grey 
lines. The shocked matter in the grey lines is hardly ob- 
served in X-rays, but the shock possibly accelerates elec- 
trons to be of order of GcV, which can be responsible 
for GHz synchrotron radio; the condition after the break- 
out is favorable for particle acceleration, as mentioned in 
Section 4.1. 
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Fig. 5. Contours of the X-ray emission measure in 
the equator (left) and polar (right) directions of 
model B2 are drawn by black lines linearly from 
zero, every 1 X 10 19 amu 2 cm -5 to the maximum 
~ 1.3 X 10 20 amu 2 cm -5 at 980 yr, every 5 X 10 17 amu 2 cm -5 
to the maximum ~ 3.8 X 10 18 amu 2 cm -5 at 1800 yr, 
and every 2 X 10 16 amu 2 cm -5 to the maximum 
~ 1.9 X 10 17 amu 2 cm -5 at 10000 yr, after explosion. 
The grey contours represent the low-level emission measure 



drawn linearly from zero, every 1 X 10 amu cm 



up to 

2 r>m " 



amr cm 



1.0 X 10 17 amu 2 cm" 5 at 980 yr, every 1 X 10 16 
up to 1.0 X 10 17 amu 2 cm -5 at 1800 yr. The horizontal and 
vertical axes show the scale in units of pc. 



5. Conclusion 



We investigate the evolution of SNRs that explode in 
the progenitors' stellar wind matter, considering possi- 
ble environments of mixed-morphology SNRs with over- 
ionized plasmas. We summarize the results; 
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• When the blast wave breaks out of the wind matter 
into the ambient interstellar medium, the shocked 
matter cools rapidly due to adiabatic expansion. 
Just after the break-out, the expanding velocity be- 
comes faster by a factor up to two, and then grad- 
ually decreases to that of the extrapolated from the 
velocity trend before the break-out. 

• Before the break-out, the shocked matter reaches 
ionization equilibrium and equipartition T c ~ J], but 
deviates from equilibrium by rarefaction after the 
break-out. Consequently, the shock-heated ejecta 
turns to be a recombining plasma, since cooling due 
to adiabatic expansion is much faster than recombi- 
nation. 

• The recombining state of the shocked ejecta lasts 
until the second reverse shock, which occurs by the 
interaction with the interstellar medium, propagates 
inward and reheats the ejecta. If the density of the 
ejecta is too low to establish ionization equilibrium, 
however, the recombining state lasts longer. 

• After the break-out in the adiabatic phase, since 
the emission measure of the shocked ejecta is much 
larger than that of the shocked ambient matter, the 
SNR in X-ray wavelengths appears much brighter 
in the reverse-shocked inner region than the blast- 
shocked outer shell. 

• When the stellar wind matter is not isotropic but 
denser in the equatorial direction due to the pro- 
genitor's rotation, the SNR in the recombining state 
looks bar-like with wings in the equatorial view and 
thin shell-like in the polar view. So that, the SNR 
would show center-filled various shapes in X-rays, 
depending on the viewing angle. On the other hand, 
the blast-shocked matter, which is very faint in X- 
rays but could be observed in radio, forms a fairly 
complete shell outside. 

• As the SNR age increases, however, the second re- 
verse shock sweeps the bar/wing structure out and 
merges into the whole ejecta eventually. Hence, the 
bar/wing structure disappears and the late-phase 
SNR would look shell-like almost independently of 
the viewing angle. 
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Midori Ozawa for meaningful discussion about Suzaku ob- 
servations of recombination radiation from SNRs. KM 
and KK were respectively supported by the Grant-in- Aid 
for Scientific Research 22540253 and 2054019, from Japan 
Society for the Promotion of Science (JSPS). 

Appendix 1. Electron temperatures of SNRs 

In the present work, we assume that shocks heat sub- 
stantially ions at the shock front, and then the elec- 
tron temperature rises through the energy transport by 
Coulomb collisions and diffusion in the post shock region. 
In Figure 6 we compare the electron temperature thus ob- 
tained (Masai 1994) by the solid lines and temperature in 



the equipartition with ions at the shock front by the bro- 
ken lines, with the observed electron temperatures (the 
circles). One can see the transfer by Coulomb collisions 
from ions to electrons is a reasonable assumption to ac- 
count for the observed data. 



i i 
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Fig. 6. Electron temperature vs. radius of SNRs which have 
no active central sources. Open and filled circles represent 
shell-like and mixed- morphology SNRs, respectively; line-con- 
nected two circles mean the data by two- temperature analysis. 
The thin and thick solid/broken lines represent the calcula- 
tions for the density p = 0.1 amu cm -3 and p = 1 amu cm -3 , 
respectively, of the ambient matter. 

References. — Kepler; Kinugasa & Tunemi 1999, G15.9+0.2; 
Reynolds et al. 2006, G27.4+0.0; Gotthelf & Vasisht 1997, 
Cygnus Loop; Miyata et al. 2007, G109. 1+1.0; Sasaki et al. 
2004, Cassiopeia A; Kawasaki et al. 2005, Tycho; Hwang & 
Gotthelf 1997, G156.2+5.7; Katsuda et al. 2009, Puppis A; 
Tamura 1995, G272.2-3.2; Harrus et al. 2001, G299.2-2.9; 
Park et al. 2007, RCW 86; Rho et al. 2002, SN 1006; 
Yamaguchi et al. 2008, Lupus Loop; Park et al. 2009, 
CTB37A; Aharonian ct al. 2008, CTB37B; Nakamura et al. 
2009, G349.7+0.2; Lazendic et al. 2005, W28; Kawasaki ct 
al. 2005, W44; Kawasaki ct al. 2005, 3C400.2; Saken et al. 
1995, Kes 27; Enoguchi et al. 2992, MSH 11-61A; Rho & Petre 
1998, 3C391; Kawasaki et al. 2005, CTB 1; Rho & Petre 1998, 
W51C; Koo et al. 1995, CTA 1; Rho & Petre 1998, W63; Rho 
& Petre 1998, HB21; Rho & Petre 1998, IC443; Yamaguchi 
et al. 2009, Kes 79; Rho & Petre 1998, HB3; Rho & Petre 
1998, G327.1-1.1; Sun et al. 1999, W49B; Ozawa et al. 2009, 
3C397; Rho & Petre 1998, MSH 11-54; Rho & Petre 1998 
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